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Binary collisions between ions and electrons in an external magnetic field are considered in second- 
order perturbation theory, starting from the unperturbed helical motion of the electrons. The 
calculations are done with the help of an improved BC which is uniformly valid for any strength of 
the magnetic field and the second-order energy and velocity transfers are treated in the interaction 
in Fourier space without specifying the interaction potential. The energy transfer is explicitly 
| calculated for a regularized and screened potential which is both of finite range and less singular than 

the Coulomb interaction at the origin and as the limiting cases involves the Debye (i.e., screened) and 
Coulomb potentials. Two particular cases are considered in detaile: (i) Ion motion parallel to the 
magnetic field with an arbitrary strength. The energy transfer involves all harmonics of the electron 
cyclotron motion, (ii) The ion arbitrary motion with respect to the strong magnetic field when 
the electron cyclotron radius is much smaller than other characteristic length scales (e.g., screening 
length, pitch of electron helix etc.). In the latter case the energy transfer receives two contributions 
which are responsible for the electron guiding center and cyclotron orbit perturbations. 
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P-f I. INTRODUCTION 

GC 

• In the presence of an external magnetic field B the problem of two charged particles cannot be solved in a closed 

form as the relative motion and the motion of the center of mass are coupled to each other. Therefore no theory 
exists for a solution of this problem that is uniformly valid for any strength of the magnetic field and the Coulomb 
force between the particles. The energy loss of ion beams and the related processes in a magnetized plasmas which 

i— i" are important in many areas of physics such as transport, heating, magnetic confinement of thermonuclear plasmas 
and astrophysics are examples ofphysical situations where this problem arises. Recent applications are the cooling 
of heavy ion beams by electrons [If, El 01 an d the energy transfer for heavy-ion inertial confinement fusion (ICF) (see, 
. e.g., Q for an overview). The classical limit of a hydrogen or Rydberg atom in a strong magntic field also falls in this 

H ' category (see, e.g., Q and references therein) but in contrast to the free-free transitions (scattering) the total energy 

i—\ , is negative there. 

Numerical calculations have been performed for binary collisions (BC) between magnetized electrons @, 0] and for 
collisions between magnetized electrons and ions [H H [l(| • In general the total energy W of the particles interacting in 
a magnetic field is conserved but the relative and center of mass energies are not conserved separetely. In addition, the 
presence of the magnetic field breaks the rotational symmetry of the system and as a consequence only the component 



of the angular momentum L parallel to the magnetic field Lu is a constant of motion. A different situation arises for 
. _ ' the BC between an electron and heavy ion. As an ion is much heavier than an electron, its uniform motion is only 
weakly perturbed by collisions with the electrons and the magnetic field. In this case Lh is not conserved but there 
j^H * 1 exists a conserved generalized energy K [ToL [TlT ] involving the energy of relative motion and a magnetic term. The 
seemingly simple problem of a charged particle interaction in a magnetic field is in fact a problem of considerable 
• • i complexity and the additional degree of freedom of the cyclotron orbital motion produces a chaotic system with two 
. £h ! degrees (or one degree for heavy ions) of freedom [l2|, [H, EH ■ 

In this paper we consider the BC between electrons and heavy ion treating the interaction (Coulomb) with the ion 
' as a perturbation to the helical motion of the magnetized electrons. This has been done previously in first order in the 
ion charge Z and for an ion at rest [l5| and in up to 0(Z 2 ) for uniformly moving heavy ion [THHy]. In Ref. [H[ three 
regimes are identified, depending on the relative size of the parameters a (the cyclotron radius), s (the distance of the 
closest approach), and S (the pitch of the helix). In earlier kinetic approaches U 0, 1 only two regimes have been 
distinguished: Fast collisions for s < a, where the Coulomb interaction is dominant and adiabatic collisions for s > a, 
where the magnetic field is important, as the electron performs many gyrations during the collision with the ion. The 
change AEi of the energy of the ion has been related to the square of the momentum transfer Ap, which has been 
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calculated up to 0(Z). This is somewhat unsatisfactory, as there is another 0(Z 2 ) contribution to AEi, in which the 
second-order momentum transfer enters linearly. Moreover, for applications in plasma physics (e.g., for calculation 
of the ion energy loss in a magnetized plasma) one calculates the angular avereged energy transfer which vanishes 
within first-order perturbation theory due to symmetry reasons and the ion energy change receives contribution only 
from higher orders Indeed, the transport phenomena, etc., are of order 0(Z 2 ) in the ion charge. 

In this paper we consider BC between ion and electrons in the presence of a magnetic field within the second order 
perturbation theory. The present paper is a continuation of our earlier study in Ref. [ll[ where the second-order 
energy transfer is calculated with the help of an improved BC treatment which is uniformly valid for any strength 
of the magnetic field and does not require the specification of the interaction potential. The paper is organized as 
follows. In Sec. |TT] starting from the exact equation of motion of two charged particles moving in a magnetic field 
we discuss some basic results of the exact BC treatment for the energy and velocity transfers as well as the energy 
conservation. In the following Sec. IIII1 we discuss the velocity and energy transfer during BC of magnetized electrons 
with ions for arbitrary magnetic fields and strengths of the electron-ion interaction potential. We assume that the 
ion mass M is much larger than the electron mass m. The equations of motion are solved in a perturbative manner 
up to the second order in Z starting from the unperturbed helical motion of the electrons in a magnetic field. Then 
in Sees. IIVI and [V] we turn to the explicit calculation of the ion second order energy transfer. For further applications 
(e.g., in cooling of ion beams) we consider the regularized and screened interaction potential which is both of finite 
range and less singular than the Coulomb interaction at the origin and as the limiting cases involves the Debye (i.e., 
screened) and Coulomb potentials. In Sec. IIVI the theory is applied to the energy transfer of an heavy ion moving 
parallel to the magnetic field B but without any restriction on B. The obtained energy transfer involves all cyclotron 
harmonics of the electron helical motion. The case of the strong magnetic field and arbitrary motion of the ion with 
respect to B is derived in Sec. [V] It is shown that the energy transfer contains two terms which are responsible for 
the electron guiding center and cyclotron motion perturbations. In Sec. lVII the results are summed up; some formulas 
for the second-order treatment are presented in the Appendices [S] and [Bj 

II. BINARY COLLISION FORMULATION. GENERAL TREATMENT 

Below we discuss the general equations of motion for two charged particles moving in a homogeneous magnetic 
field and the remaining conservation laws. From the velocity transfer we then proceed the energy transfer of particles 
during binary collision process. As shown in Ref. [Tlj |. t ne present treatment becomes more transparent in Fourier 
space. 

A. Relative motion and conservation laws 

We consider two point charges with masses mi, mi and charges q%e, qi e, respectively, moving in a homogeneous 
magnetic field B = Eh. We assume that the particles interact with the potential qiq2^ 2 U(r) with fi 2 — e 2 /4ireo, 
where £o is the permittivity of the vacuum and r = ri — r2 is the relative coordinate of colliding particles. For 
charged particles the function U(r) can be expressed, for instance, by the Coulomb potential, C/c( r ) = V r - I n plasma 
applications the infinite range of this potential is modified by the screening, e.g. UpM = e~ r / x /r with a screening 
length A which can be chosen as the Debye screening length Ad, see, for example [17| • The quantum uncertainty 
principle prevents particles from falling into the center of these potentials. In a classical picture this can be achieved 
by regularization at the origin U R (r) = (l - e- r l x ) e- r l x /r, see for example [18|, [19|, where A is a parameter, which 
may be related to the de Broglie wavelength. 

In the presence of an external magnetic field, the Lagrangian and the corresponding equations of particles motion 
cannot, in general, be separated into parts describing the relative motion and the motion of the center of mass with 
velocities v, V cm and coordinates r, R cm , respectively (see, e.g., [6|). Introducing the reduced mass = 1/ mt+l/mz 
the equations of motion are 

V(t) + fi 4 [V(t) X b] = -As [V cm (t) X b] + (r(i)) , (1) 

M 

V CBl (t) - Oi [V cro (i) x b] = -0 2 [v(i) x b] , (2) 

where qxq^^Y (r(t)) (F = —dU/dr) is the force exerted by the particle 2 on the particle 1. (The force which acts 
on the particle 2 is — qiq2^ 2 F (r(t))). The frequencies fii, f^, ^3 and SI4 are expressed in terms of the cyclotron 
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frequencies u> c i = \q\\eBjm,\ and w c2 — |<7 2 |eJ3/m 2 of the particles 1 and 2, respectively, 



mi?io; cl + m 2 <; 2 w c2 „ ^2^ c2 - ftWcl ... 

S2l = ■ , il 2 = ■ [I, (3) 

771 1 + m 2 mi + m 2 



S2 3 = ? 2 a> C 2 - ftWd, S2 4 = -A* 1 ■ ( 4 ) 

777i 7772 



Here <r„ = /q v with i/ = 1,2. From Eqs. ((TJ and ([2]) follows the conservation of total energy 

w = i + ^ + H<l2^U(r) = const, (5) 

but the relative and center of mass energies are not conserved separately. 

The coupled, nonlinear differential equations (JT|) and @ completely describe the motion of the particles. They have 
to be integrated numerically for a complete set of the initial conditions for solving the scattering problem. In the case 
of interaction of heavy ions (m 2 = M, g 2 = Z) with electrons (mi = m, q\ = —1), i.e. M 3> m, the equations of 
motion can be further simplified, since (J, — > m, fix, fi 2 — * and ^3,^4 — > w c (see Eqs. (J3J) and Q), where w c = eB/m 
is the cyclotron frequency of electrons. Equation © leads to V cm — > Vj = const, where Vj is the heavy ion velocity, 
and Eq. (JTJ) turns into 

v(t) + w c [v(i) x b] = -cu c [ Vi x b] - -£-F (r(t)) . (6) 

m 

With the help of the equation of motion ([B]) it can be easily proven that the quantity 

2 

is a constant of motion. In contrast to the unmagnetized case, it thus follows that the relative energy transfer during 
ion-electron collision is proportional to Srj_Vij_, where 5r±_ and Ujj_ are the perpendicular components of the change 
of relative position and the ion velocity. 

B. Energy loss and velocity transfer 

The rate at which the energy of an ion in a collision with an electron at time t changes is given by 

dEi(t) 



dt 



^vrF(r(f)), (8) 



as Z^ 2 F (r) is the force exerted by the electron on the ion. Integration with respect to time yields the energy transfer 
itself 

5E l {t) = Zf f Vj • F (r(r)) dr , (9) 

J — OO 

which after completion of the collision becomes 

/oo 
Vj • F (r(r)) dr. (10) 
-oo 

Here one has to integrate the equations of motion for the relative trajectories. The limit M 3> m leading to Eq. © 
implies that the change in the ion energy is calculated under the assumption of a constant ion velocity. Alternatively 
this energy transfer can be expressed by the velocity transferred to the electrons during the collision. For that purpose 
we substitute Vj = v e (t) — v(t) into Eq. (Tl0|) and split the electron velocity into two terms v e (t) = v e o(i) + Sv(t), 
where v e o(t) describes the helical motion in the magnetic field 



V e0 + LU C [v e0 X b] = 



(11) 
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and 8v(t) the velocity transfer (we assume that 8v(t) — * at t — > — oo) due to the collision with the ion 



8v(t) + lo c [<Jv(t) x b] = -^-F (r(t)) . 



This yields 



/ v e0 (r) -F(r(r))dr+ f 5v(r ) • F (r(r)) dr+U (r(i)) 

J — oo J —oo 



(12) 



(13) 



In Eq. (fT3")l we take into account that at t — > -co, r(f) — » oo and U(r(t)) — ► 0. The time integrals in Eq. (T5j) can be 
done with the help of the derivative of the scalar product v e0 (£) • 5v(t). Using the equations of motion ([IT]) and (fT5|) 
we obtain 



[veo(t) ' «v(t)] 



v e0 («)-F(r(*)), 



which yields 



Similarly from Eq. (fT2"|) 



which yields 



Thus 



1 f v e0 (r) • F (r(r)) dr = -mv e0 (t) • Sv(t). 

J —OO 



Z<j? [ 5w{t) -F(r(r))dT 

«/ — OO 



-?[Mt)] 2 - 



= ^ 2 C7 (r(t)) - mveo(i) • Sv(t) 
The last two terms in this equation represent the change in the electron energy due to the collision 

5E e {t) = y {[veo(i) + 5v(t)f - v e0 (t)} = mveo(i) • 8v(t) + y [Sv(t)f . 
This shows energy conservation 

8Ei(t) + 8E e {t) - Z</?U{r{t)) = 0. 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



As U(r(t — > oo)) = the energy change of the ion can also be calculated from the velocity transfer Av = 8v(t — > oo) 
with the help of 



AE, = -<>•£•. ( / - - oc i = -m [ Veo(t) • Av + ^Av z j . 



(21) 



this method has been adopted in jTy] . In this approach the potential U(r) has to be specified at an early stage. In 
Sec. IIIII we will show that Eq. [)10p allows for a more general formulation in which the cut-off at large distances and 
the regularization at small distances can be treated easily. 

Until now we have considered the energy transfer of an ion. In addition this energy transfer, AEi, can be expressed 
by the change of the relative energy, AE r , and the electron momentum transfer Ap = mAv. Because we are dealing 
with heavy ion with Vj = const, the relative velocity transfer is the same as for the electrons, i.e. 8v(i). Since the 
unperturbed relative velocity is vo(t) = v e o(t) — Vj (see Eq. (p7)l below) we can establish a simple relation between 
energy transfers AEi and AE r given by 



AEi = -AE r - Vj • Ap. 



(22) 



It is clear that the relative energy and momentum transfers depend only on the relative quantities and the ion velocity 
Vj is not explicitly involved in AE r and Ap. Thus, having the ion energy transfer the other quantities can be easily 
extracted from Eq. 
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III. PERTURB ATIVE TREATMENT. GENERAL THEORY 
A. Trajectory correction 



In this section we seek an approximate solution of Eq. (fl2| in which the interaction force between the ion and 
electrons is considered as a perturbation. Thus we have to look for the solution of Eq. (fT2)l for the variables r and v 
in a perturbative manner 

r(f)=r (t)+ri(t) + r a (i)..., v(t) = v (t) + vi(*) + v 2 (*)..., (23) 

where i"o(i),Vo(i) are the unperturbed ion-electron relative coordinate and velocity, respectively, r n (t),v n (t) oc 
Z"F„_i (n = 1,2,...) are the nth order perturbations of r(t) and v(t), which are proportional to Z n . F n (t) is 
the nth order correction to the ion-electron interaction force. Using the expansion (|23|) for the nth order corrections 
F„ we obtain 

F (r(t)) = F (r (t)) + F x (r (t), ri(t)) + (24) 

where 

F (r (t)) = F (r (t)) = -i / dk£/(k)kexp [tk • r (i)] , (25) 

Fi (r (t), r x (t)) = (n(t) • |-) F(r) = / dk£7(k)k [k • n(t)] exp [ik • r (t)] . (26) 

V cry r=ro(t) J 

In Eqs. (]25fl and l|26[). we have introduced the ion-electron interaction potential U(r) through F(r) = —dU(r)/dr and 
the force corrections have been written using a Fourier transformation in space. 

We start with the zero-order unperturbed helical motion of the electrons. From Eq. Ijlip we obtain 

vo(i) = v r + v e ± {ucos(uj c t) + [b x u] sin(w c i)} , (27) 

ro(t) = Ro + v r i + a {usin(w c i) — [b x u] cos(w c i)} , (28) 

where u = (cos <p, sin tp) is the unit vector perpendicular to the magnetic field, v e » and v e ± (with v e ± > 0) are the 
electron unperturbed velocity components parallel and perpendicular to b, respectively, v r = u e ||b — Vj is the relative 
velocity of the electron guiding center, and a = v e j_/u c is the cyclotron radius. It should be noted that in Eqs. (|27p 
and (|28[) . the variables u and Ro are independent and are defined by the initial conditions. 
The equation for the first-order velocity correction is given by 

vi (i) + uj c [vi(t) x b] = -^F (r (t)) (29) 

TO 



with the solutions 



vi (*) = ^-{-bV\\ (t) + Re [b (b • V ± (*)) - V± (t) + * [b x V± (t)]] } , (30) 



77) 



ri(t) = — {-bP)i(t) + Re [b (b • P ± (t)) - P ± (t) + i [b x P x (t)]]} , (31) 



f 

777 

where we have introduced the following abbreviations 



V\\ (t) = ( drb- F (r (r)) , Vx(i) - e*"»* / dr e -^ T F (r (r)) , (32) 

J — CXD J — OO 

Py (*)=!/+/' dr(i-r)b.Fo(ro(r)), P±(t) - — / dr fe^^ - ll F (r (r)) (33) 
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with v — — iVj|(t)| _ and have assumed that all corrections vanish at t — > — oo. For instance, in the unscreened 
Coulomb case, the interaction force Fo must behave as Fo (ro(f)) — > l/t 2 for |t| — > oo. Thus from Eq. (f3"2")l at t — > oo 
we obtain Vji(t) — * Voll = const and Vj_(i) — > e l " c *Voj_, where Voj_ = const. The quantities Voii and Voj_ give the first 
order velocity correction in Eq. (|30|l after an electron-ion collision. In this limit, we find for the first order trajectory 
correction from Eqs. d32J) and $33J) P\\(t) = V \\t + P \\, P±(f) = -i (V _l/u; c ) e i£ "<=* + P j_, where 

/oo • /-oo 

drrb • F (r (r) ) , P 0± = — / drF (r (r) ) . (34) 
-oo w c J_oo 

For the Coulomb interaction v — h-v r /v^ and v = for any screened interaction potential. Note that for the Coulomb 
interaction the second term in the first relation of Eq. (|34p tends to infinity (see, e.g., [HI). However the contribution 
of this term to the ion energy change vanishes after averaging over impact parameters. 

Substituting Eqs. I|25p and (|28p into Eq. (|33p . and using the Fourier series of the exponential function e lzs ™ e [20| . 
we obtain for an arbitrary interaction potential 

flift) =i / dkC/(k)(k-b)e ik - Ro V e mip J n (k ± a) — — (35) 

" J „^oo '(^(^-ZO) 2 ' 

Px(t)=i / dkC/(k)ke ikRo ^ e^J n (fc ±a ) (36) 
y n f^ oo (Cn (k) - «0) (U-i (k) - «0) 

where J„ are the Bessel functions of the nth order. Here (k) = ncu c + k • v r , tp = tp — 9, tan 8 = k y /k x and k± is the 
component of k transverse to the magnetic field. The quantities Vn(t) and Vj_(t) are obtained directly from Eqs. |35|) 
and ([36]) through the relations Vj| (t) = iVi) and V±(t) = P±(t). 

It should be noted that Eqs. p0|) and (f3Tj) give formal but exact solutions for the velocity and trajectory corrections 
to the unperturbed quantities in Eqs. (|27|) and (|28|) if the first order force Fo (ro) in Eqs. ([32| and (|33|) is replaced 
by the exact one, F (r) with exact relative coordinate r. Substituting Eqs. (f3"0"P and (pJTj) with exact force F (r) into 
Eq. (|2ip we obtain an exact relation for the ion energy transfer 

AE i = Z^ 2 {V||« e || + v e± [(u ■ V c ) - V. • [u x b]]} (37) 



Z 2 ^ 
+ — — 

2m 

Here 



{2 (V c • [b x V,]) - Vf - [b x V c ] 2 - [b x V s 



2 



V„ = dr h • F (r(r)) , { £ = £ drF ( r(r )) { £ ^ (38) 
are constants. The exact energy transfer (|37p is now expressed by only the relative coordinate r(t). 

B. First and second order energy transfers 

The total energy change of the ion during an ion-electron collision is given by Eqs. (fTUp . Insertion of Eq. (f2"4"P into 
the general expression (fTUp yields 

AEi = AE^ + AE^ + (39) 

where 

/oo ^00 
cftv; • F (r (i)) , A£,f } - Z^ 2 / rftv, ■ Fi (r (t), n(t)) (40) 

are the first- and second order energy transfer, respectively. 
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1. First order energy transfer 

The first-order energy transfer can be obtained by substituting Eqs. (|25|) and (|28]) into the first one of Eqs. ((40)) . 
This yields 

/oo 
dkU{k) (k-vOe*- 110 ^ e mv> J n (fcj.a)<J (C„(k)) . (41) 
n— — oo 

We now introduce the variable s = R _/ which is the component of Ro perpendicular to the relative velocity vector 
vy. From Eqs. (|27|) and ([28]) we can see that s is the distance of closest approach for the guiding center of the electron 
helical motion. The energy loss is now given by the average of AEi with respect to the initial phase of the electrons 
ip and the azimuthal angle of s. For spherically symmetric interaction potentials (U(r) = U(r) and U(k) = U(k)) the 
first order energy transfer gives no contribution due to symmetry and the ion energy change receives a contribution 
only from higher orders. In fact, Eq. (|41[) for the averaged first order energy change gives 

(A£ 4 (1) ) = -2mZf J dkU(k) (k • v,) J (k,s) J (k ± a)S (k • v r ) , (42) 
where k 2 = k 2 — (k • n r ) 2 and n r = v r /v r . As the integrand is an odd function of k we have (AE^ 1%> ) = 0. 



2. Second order energy transfer 
Inserting Eqs. (US]), fl25J), ([5TJ). l[55]) and flSgJ) into the second equation of Eqs. (gUJ) one obtains 

A R (2) (Ro, <p) = [ dkdk'U{k)U(k') (k • vO e < k+k ')- R o (43) 

TO J 

J- e i ^+ i ^'j n (k ± a)J m (k' x a)S(<: n (^) + Cm(k')) G m (k, k'), 



+ OC 

X 

n:m=— oo 



where ij;' = cp — 8 f , and 



ggi (k,k') g2 (k,k')-i g3 (k,k') 

mi ' j ~ (C m (k') - *o) 2 + (C m (k') - «o) (C™-i(k') - io) (44J 

ff2 (k,k')+i<?3 (k,k') 



(C m (k')-i0) (C m+ i(k')-i0) 



with 

ffi(k,k') = -(k-b)(k'-b), g2 (k,k') = (k-b)(k'-b)-(k-kO, (45) 
53 (k,k') = k-[k'xb]. 

(2) 

Next, for the practical applications, AE\ ' is averaged with respect to the initial phase of electrons <p and the 

(2) 

azimuthal angle $ s of the impact parameter s. After averaging the energy transfer AE\ with respect to ip the 
remaining part will depend on 5 ((k + k') • v r ), i.e. the component of k + k' along the relative velocity n r . Thus this 

(5-function enforces k+k' to lie in the plane transverse to n r so that e*( k+k )' R °(5 ((k + k') • vy) = e lC *' s 5 ((k + k') • v r ), 

where Q = k^ + k'^' and k^' is the component of k transverse to n r . The result of the angular averaging reads 

(AEW)^. = [ dkdk'U{k)U{k') (k ■ Vi ) J (Qs) 5 ((k + k') • v r ) (46) 

TO J 

OO 

]T (-l)"e m ( 9 - 9 ')j„(fc ia )J n (^a)G„(k,k'). 



x 

n— — oo 



This series representation of the second-order energy transfer is valid for any strength of the magnetic field. Besides 
the direction b of the magnetic field and the direction n r = v r /v r of the relative velocity is singled out in the argument 
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of the 5-function and summand of the n-summation. This prevents a closed evaluation of the energy transfer (|46|) . 
However the limiting case of an ion motion parallel to a magnetic field of arbitrary strength and the case of an 
arbitrary motion in a strong magnetic field can be treated in a straightforward manner. 

The calculation of the angular-averaged relative energy and momentum transfers is performed by the similar method 
as for deriving Eq. (|46[) . Note, however, that due to the symmetry reason the transverse part of the momentum transfer 

vanishes, v i± • (Ap ± ) = 0, and Eq. J22J) is reduced to (AE^) = -(AEi 2) ) - vq(Ap^). 

It is also usefull to integrate the <p, $ s -averaged ion energy change, (AE> ) ip ^ s , over the impact parameters s in 



the full 2D space. Thus we can introduce a generalize cross section 



I El 



through the relation 



a(v r , Vl ) = 2ir (°°{AEf = - ^ ^ / dk|£/(k)| 2 (k • v 
Jo 2mv r J 



2mv r 

2 



(47) 



J2 Jn(k±a) k\8' (Cn(k)) + ^-[6 (C„+i(k)) - S (C„-i(k))] f , 



n— — oo 



where S'(x) defines the derivative of the 5- function with respect to the argument. Note that for the Coulomb interaction 
U(k) = Uc(k), the full 2D integration over the s-space results in a logarithmic divergence of the k-integration in 
Eq. (|47|) . This will be explicitly shown in the next section. To cure this, we should introduce the cutoff parameters 
k mm and fe max , see [ll[ for details. 



IV. ION PARALLEL MOTION AND ARBITRARY MAGNETIC FIELD 



The averaged energy transfer, Eq. (I46p . can be evaluated without further approximation for an ion motion parallel 
to the magnetic field and assuming axially symmetric interaction potential, U(k) — ?7(|fe|| |, fcj_). In this case the 
averaged energy transfer can be represented as the sum of all cyclotron harmonics. For the parallel motion of the ion, 
Vjj_ = and v r = (z; e || — v.^ b = iV||b, where i^m and Vjj_ are the components of ion velocity along and transverse to 
the magnetic field, respectively. In general case setting Vjj_ = 0, we have from Eq. (|46[) 



T7iZ 2 ^V, 



e v-n f 

f—p- / dkdk'l7(k)l7(k / )fc||Jo(|k / x + kj.|a)5(ifcf|+A!||) 
x £ (-l)"e m ( e - e ')j„(fc x a)J„(fc ±a )G„(k,k'). 



rn | ■ 

oo 



(48) 



Introducing cylindrical coordinates along b and using the addition theorem for the Bessel functions Jo (|k^_ + kjjs) 
[20T | the energy transfer (|4"8|) finally reads 



(2), i&fVi 



7 ' p v i\\ 2 { / \ ® t 

—3- 2^ n i 3C/„ (A;,, , a, s) + fey frj-Un (fe|| , a, 

r\\ 11=1 "II 



+ — [V n (fc|| , s, a) - V n +i (fey ,s,a)+ V n (k\\ , a, s) - V n+ i (fe|| ,a,s)] 
where 6 — \v r u\/uj c is the pitch of the electron helix, divided by 2ir, and 



kn —nj 8 



(49) 



U n (k\\,a,s) = 
V n (k\\,a,s) = 



(2nf 



(27T) 



/•oo 

/ U (fe||,fej_) J n (kxa) J n (k±s) kj_dkx 
Jo 

/>OC 

• / U (fe||, fcj_) J n (k±a) J n -i (fcj_s) kj_dk± 
Jo 



(50) 



(51) 



We recall that for the ion parallel motion the last term in Eq. ([7]) vanishes, i.e. the relative energy is conserved 
(2) 

and (AEf. ) = 0. Therefore, the relation between ion energy and relative momentum transfers is simplified to 
(AE& ) = -v i]l (Ap\V). 
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Now we specify the electron-ion interaction. In the following we consider the regularized screened potential U (r) 
U-r(t) introduced in Sec. Ill Al with 



(2tt) 2 \kl + K 2 kl+x< 



(52) 



where k 2 — fcjj + A 2 , x 2 = k 2 + d 2 and d 1 = A 1 + A . Carrying out the calculation of the fc^-integrals in 
Eqs. (|5T)|) and (f5"Tj) (see, e.g., [2(|) with the potential (p)2")) and substituting into Eq. I|49p for the regularized screened 
interaction we obtain 



(AE, 



(2)\ 



AZ 2 </> i v 



^ n 2 < 3 [u„ {n n a, n n s) - u n (x n a, Xns) 

r\\ n=l { 



2n 2 



—~T n {Una, K n s) -T n {Xnd, XnS) 

n n An 



(53) 



Here 



+ 2S 2 [n 2 n Q n (k„o, k„s) + xlQn (Xn«, Xn,s) - 2k„x„-D„ {K n a, n n s; x n a, Xns)] 



n 2 1 n n 1 
i-.^ = i = i 

S 2 A 2 ' (S 2 d 2 ' 



(54) 



Qn(x,y) = I n (£)K n (ri) 



(55) 



with £ = min(:r, y), 77 = max(:c, y), and the modified Bessel functions /„ and X n , 



D„(a;,y;X,y) = — [s n -i (x,y) s„_i (-X",Y) - s„ (2,2/) s„ (X,Y) 
4n 

+ s n -i (y, z) s„_i (Y, X) - s n (y, x) s n (Y, X)} , 



(56) 



!I n (x) K n+1 (y) y>x 
-In+i (y) K n (x) y < x . (57) 

\ [l n (x) K n+1 (x) - I n+1 (x) K n (x)} y = x 

For a study of the convergence of the series in Eq. (|53|) we note that in all terms the modified Bessel functions /„ carry 
the smaller argument oc min(a, s), while the K n depend on max(a, s). At large harmonic numbers n both the indices 
and the arguments of these functions are large, and they behave as I n (nt;), K n {nrj). Therefore the case s — a is most 
critical for the convergence of (|53[) . This is intuitively clear as the gyrating electron hits the ion on such a trajectory. 
This should not matter for the potential (|52p. which has been regularized near the origin for exactly that purpose. 
Since at large n — > 00 and at s — a the summand in Eq. (|53p involves the functions I n (nx) and K n (nx) and their 
derivatives the further analysis can be done using the uniform asymptotic expansions of the modified Bessel functions 
[20I I2H ]. Insertion of the expansions shows indeed that the nth member of the series are of the order 0(n -4 ), so the 
series converges even for s — a. On the other hand, the energy transfer for the unregularized potentials Uq and {7 D 
diverges for s = a. 

For the screened but unregularized potential, i.e. in the limit A — > 0, all functions u n , T n , Q n and D n involving Xn 
in their arguments tend to zero. There remains 

< A 4i } > - -^kr £ " 2 [3m - (K - a ' KnS) (58) 

r|| n=l 

2n 2 2 
+ 1 n 7 i u n K n s ) T n (n n a, K n s) + 2(n n 6) Q n (n n a, k„s) . 
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Now insertion of the uniform expansions of the modified Bessel functions shows that the members of this series are 
independent of n for n — ► oo and s = a, hence the series diverges. In the Coulomb case, i.e. for A — ► oo, K n = n/8, 
the resulting series are geometric and can be summed in closed form 



'I " 2m5 2 vf l{ + + v 2 ) sinh*(£, n) 



(59) 



x [2 + S 1 (£,r y )]e-*^ + 



sinh*(^,r?) 

with f = min(f , f), 57= max(§, §), *(£,??) = - ¥>(r?) and 



V (0 = VeTT In [ V ^t I+1 ] , (60) 



5(l + e 2 ) 3/2 5(1 + 7y 2 ) 3 / 2 



6£ 2 (l + ?7 2 ) 3 / 2 6t7 2 (1+£ 2 ) 3 / 2 

(l + ?7 2 ) 3 / 2 (1 + C 2 ) 3/2 

2r/ 2 (l+C 2 ) 1/2 _ 2£ 2 (l + r? 2 ) 1 /2' 

(1+^2)3/2 (1 + ^2)3/2 



3 U 2 



(61) 



For the limit \s — a| — ► the Taylor expansion of the function 77) is used. This yields 



( A 4i > " 



1 



1 



m ^' y r|| ^(^ + r; 2 ) 3 / 2 |s — a| 



(62) 



which exhibits the divergence at s = a. Note that in Eq. (|62|) 77 = a/<5 = v e ±/|iV||| and does not depend on the 
strength of the magnetic field. 

For later purposes we also note the limits of Eq. (|58p for a small electron transversal velocity with a <C 6, 



2v, 



rnv . 



1 + (/C]A KlS, K2S) 



(63) 



with 



2 



XKi (A) V^iPO 



1 

C 2 " 



Ai^ (A) 



^ (A) 



(64) 



Because of the symmetry of Eq. (|58jl in respect to its arguments the limit of a small impact parameter s <C 5, is given 
by Eqs. (|6"3")l and ([M]) with the roles of a and s interchanged. Similarly from Eq. ([53")) in the case of the regularized 
potential and for vanishing cyclotron radius (v e ± = 0) we obtain 



2i-i 



j-KKiaJifiM-OtiaJ/iriOci*)] 2 - 



(65) 



Comparing the first term in Eq. (|63p with Eq. I|65p we conclude that the additional modified Bessel function in Eq. 
with the argument xi s guarantees the convergence of the energy transfer at small impact parameter s. 

For the practical applications and for general interaction potential in Appendix[X]we also performe the s-integration 
of the second order energy transfer, Eq. (|49| . 



V. ARBITRARY ION MOTION IN A STRONG MAGNETIC FIELD 



After the discussion of the energy loss for an ion moving parallel to a magnetic field of arbitrary strength we return to 
the general case, where the ion velocity has a component transverse to the field. As mentioned above the integrations 
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and summation in Eq. (|46[) cannot be done in this case unless other simplifications are made. In the following we 
consider strong magnetic fields and in Eq. (|46[) keep only the terms with n = 0, ±1 and expand the remaining Bessel 
functions Jo(k±a) ~ 1 — (k±a) 2 /4:, Jx{k±a) ~ k±a/2 with respect to the cyclotron radius of electrons. Note that this 
approximation can be alternatively formulated as a smallncss of the electron transverse velocity v e ± . We obtain two 

(2) (2) 

contributions to the energy transfer, (A£7- )i which is independent of the cyclotron radius and (AE^ }n which is 
proportional to a 2 . We split (AE^}i according to 

(AeI% = {Ae(% 1 + {Ae(% 2 , (66) 

where 

(AE^)n = (2 ^ 2r? f^ / dkdk'U(k)U(k') (k ■ v,) (k • b)(k' • b) Jo (Qs) (67) 
x <5((k + k') . Vr )<J'(k-v r ), 



(A^ (2) )i 2 - /dkdk'C/(k)C/(k')(k- Vi ) J (Qs)<5((k + k')-v r ) 

muj c J 



(68) 



ff 3 (k,k') 



1 



1 



k • v r k • v r — uj c 



7rg2 (k, k') S (k • v. r — u) c 



As in Sec. IIVI we will assume axially symmetric potentials U(k) = U(\k\\\, k±). The integration is done by using 



cylindrical coordinates oriented along n r = v r /w r , i.e. any vector C will be represented as C = Cijj n r + C^' . For 
the Bessel functions we use the addition theorem [201 ] . The angular integrals are trivial as they involve powers of 
trigonometric functions. Note that the contribution of the term proportional to the function ^(k, k') in Eq. 



( T"l f( 7'') 

vanishes due to the antisymmetrical behavior of this function with respect to the azymuthal angles of k_^_ and kj; . 
Then for the energy transfers after straightforward calculations we obtain 

Z 2 ^v 2 ± 



(A^ (2) )n 



inv'. 



( v e\\ - v 'i ) T 12 ( s ) + v e \\v r \\Toi (s) T 03 (s) 



(69) 



(A4 2) >i 2 



Z 



2J.4 



v i± [ 2< 7 2 ( v e\\Vi\\ - v 2 ) 7^(g, s) - v e \\v r \\T 01 (q, s)T 03 (q, s)] 



(70) 



+ 



2v^(v e \\Vi\\ - V 2 ) - uf X ( U e|| ~ V i) T 12^ s) 



where q = 5 1 = uj c /v r , T„ M (s) = T vll (0, s) and 



T 

- 1 V^L 



(27T) 2 



U(q,k±)J v {k±s) k^_dk_L. 



(71) 



We now specify the interaction potential and explicitly calculate the functions ([7T|) . From Eqs. ([52"]) and ([7T|) we 
obtain in the regularized and screened case 



T u(^ s ) = K i K i ( K i s ) - XiKi (Xi s ) ■ 
^oafes) = xlKo (Xis) - n\K ( Kl s) . 
T*(q,s) = K (tiis)-K { X is) 



(72) 



with k\ = q 2 + A 2 and \i 



= q 2 + dr 2 . The functions (s), (s) and (s) are easily obtained from Eq. (J72J) 
setting there q = 0, i.e. n\ = 1/A and xi = 1/^- We investigate the asymptotic behavior of the functions in Eq. (f?2"]) 
in the limit s — ► 0. As iio (z) ~ ln(l/z) and ifi (z) ~ 1/z the divergence is not worse than logarithmic and will cause 
no harm when integrating over the impact parameter s. This is not so in the case of the screened potential. Insertion 
of Eq. ((72j) at A -> (i.e. xi -> oo) into Eq. {69]) yields 



(AS| 2) )I1 = 



— ) {^\\ ~ v i) [P K ^P)f - MM [P K o(p)} 2 } 
s J mv° l 11 j 



(73) 
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with p = s/A, which behaves like s~ 2 for small impact parameters. On the other hand all functions in Eq. (|72[) vanish 
exponentially for large impact parameters because of the finite range of the potentials C/ R and U-q . We obtain the 
Coulomb case by taking the limit A — > oo in Eq. (|73| . This yields 



(A^ (2) )n 



Z f\* „,2 



a 



K 2 || " «?) 



(74) 



which is precisely the energy transfer for tight helices |16l |. 

We turn now to the next term {AE^)\2 in Eq. (|66|) . i.e. Eq. (|70|) with functions from Eq. (|72|) . For the screened 
potential Uu the insertion of Eq. (f?2")) at A — > yields 



(A^ 2) }i 2 



(zp\' l 



- 



(f e || ^|| - vf) + V e \\V r \\ 



[piK (pi)f 



(75) 



+ 



(2) 

with pi = K\s. For parallel ion motion, Vi± — > 0, the first term vanishes {AE\ )n =0 and there remains a contribution 

to (AEj }i2, which as expected is equal to the leading term of the corresponding expansion of Eq. (|63|) in orders of 
a/5. 

(2) 

Taking now the limit A — > oo for the unscreened Coulomb interaction we see that \AE> )u vanishes exponentially 
for s > S. Hence the energy transfer is given by the tight helix term Eq. (|74|) . For the case s < 5 and hence p\ ~ 
we obtain 



(AE™ > n + <A25< a > ) I2 



2v, • v r 



(76) 



This is just the stretched helix case considered in Ref. [16j . Equation ([76)) is the same as the second-order energy 
transfer in a field- free case (see, e.g., Refs. [11]]) but the full relative velocity u r = v e — Vj is replaced here by the 
relative velocity of the electron guiding center v r . 

These results are obtained in the limit a — > where the electrons move along their guiding center trajectories. 
Moreover, for oj c — > oo also the pitch 5 — > and these trajectories are rectilinear along the lines of the magnetic 

(2) 

field and the energy transfer is given by (AE^ )n- For a finite cj c corresponding to a finite pitch the contribution 

(2) 

\AE> )i2 describes the perturbation of the guiding center trajectory. 

The quadratic term (AE^)u ~ a 2 accounts for the finite cyclotron motion of the electrons. In general this term 
is obtained from Eq. (|4*6")l and reads 



<A£f>>n = 



4m 



; J dkdk'U(k)U(k') (k • v 4 ) J (Qs) 6 ((k + k') • v P ) 



(77) 



{(Jfei + k'l) G (k,k') + k ± k' ± [e l ( e - e ') Gl (k,k') + e-*( e - fl ')G_ 1 (k J k')] } 



Here Go(k, k') and G±i(k, k') are given by Eq. (|4"4")l . Using the same techniques as before the straightforward 
calculation yields 



<A£f>) II= ^ 2a 



1 



(Qo + Qi + 0,2 



\ s 2 J 2mv 2 

Here we restrict ourselves to give the result for the screened potential Uu 



(78) 



So = 



v e\\Vr\\Vjl. 
V 4 

v e\\Vr\\ 



V ^ r H \r?K InW 2 X ^ 



U„||U. 



Pill 



3 V ? 



2w 2 



4vl 



(79) 



[p #1 (p)J + L_i ^ Jf 2 
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Qi = 



+ 




[p\Ko(pi)\ 



(80) 













Mi. 


,,,2 





[p\K2(n)Y 

4g 4 ^|| u tl 



Pi«b(pi)ifi(pi), 



U P ||U r 



V.\\Vi 



2v 2 



8q 2 



[pIkmY 



2v 



1 - 



Avl 



Aq 2 



2(, 



'e\\Vi\\ 



[pIkm]' 



(81) 




[pIK 2 {P2)\ 



Here k 2 = A/8 2 + 1/A 2 , p 2 = k^s. The cyclotron motion and the drift of the guiding center of the electron are coupled 
to each other. Therefore the perturbation of the cyclotron motion causes an additional perturbation of the guiding 
center motion. This effect is given by the first term Qq in Eq. (|78[) which depends on magnetic field through the 
cyclotron radius a in the prefactor, while the arguments of the modified Bessel functions in Eq. (|79p do not depend 
on magnetic field. In the other terms Q\ and Q2 the arguments of the Bessel functions p\ and pi correspond to the 
first and second cyclotron harmonic perturbations, respectively. 
For Vi± = we have Qo = and 



pUk^-kKp,)] 



6 



p\K 2 { Pl ) 



:P 5 1 K (p 1 )K 1 ( Pl ) 



(82) 



Q 2 = [plK 2 {p 2 )} 2 . (83) 

iV || 

With the help of the recursion relations of the modified Bessel functions it is easy to see that the resulting energy 
transfer {AE^^n agrees with the corresponding 0(a 2 /(5 2 )-term of the energy transfer, Eq. (f63|) . where the limit of 
parallel ion motion vn_ —> was taken before the limit a -C S. 



VI. DISCUSSION AND CONCLUSION 



In this paper, we have presented a detailed theoretical investigation of the energy transfer of a uniformly moving 
heavy ion due to the binary collision (BC) with the magnetized electrons. The BC energy transfer can only be 
evaluated explicitly in closed form in the limiting cases of a vanishing and an infinitely strong magnetic field. The BC 
treatment developed here is valid for arbitrary strengths of the magnetic field and arbitrary shapes of the interaction 
potential up to second order in the interaction strength. The purpose of this work was to investigate the ion energy 
transfer for finite magnetic fields which is explicitly calculated for a regularized and screened potential which is both of 
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finite range and less singular than the Coulomb interaction at the origin and as the limiting cases involves the Debye 
(i.e., screened) and Coulomb potentials. Two particular cases have been considered in detail: (i) Ion motion parallel 
to the magnetic field with an arbitrary strength. The energy transfer involves all harmonics of the electron cyclotron 
motion, (ii) The ion arbitrary motion with respect to the strong magnetic field when the electron cyclotron radius 
is much smaller than other characteristic length scales (e.g., screening length, pitch of electron helix etc.). We show 
that in the latter case the energy transfer receives two contributions which are responsible for the electron guiding 
center and cyclotron orbit perturbations. 

We would like to mention that our current results leave still some questions open. It is clear that for the validity of 
the second-order perturbation BC theory developed here more critical are the domains of the small relative velocities 
v r and/or impact parameters s. Moreover, for the binary electron-ion collisions in a magnetic field as given by the 
equation of motion ^ there are less integrals of motion than degrees of freedom which indicates the possibility of the 
chaotic dynamics in the system (l2l . [l3l Il4| . This immediately raises the question whether a perturbative treatment 
as proposed in this paper can be applied at all. We will address this issue in the fortcoming studies by showing some 
examples for the energy transfer obtained from a numerical solution of the equation of motion ([6]) and by making 
some comparison with the perturbative treatment. This topic is presently under investigation and will be published 
elsewhere. 
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APPENDIX A: INTEGRATED ENERGY TRANSFER FOR PARALLEL ION MOTION 

The integration of the energy transfer Eq. (I49| with Eqs. (f50|) and (f5Tj) with respect to the impact parameter s is 
facilitated by using the following relations for the Bessel functions 

/>oo />oo 

/ J n (k±s) J n (k' ± s) sds = / J n _i (k_Ls) J n -i {k' x s) sds = — S (k' ± - k±) . (Al) 
./.. Jo k: 

Using this relation we integrate the energy transfer Eq. (|49|) with respect to the impact parameter. This yields 
Z* 00 f 7 J* 2 \ 4 00 f r) 

jf (AE^)sds = ^■En a { 3 *»(*B.°) + *ya3bjj-*n(*y' ) ( A2 ) 
s 2 1 

+ 2^[*»-i(*l|.o)-*«+i(*ll»o)]] 



where 



$ n (fc||,a) = / U n (fey, a, s) sds = / U 2 (k h kj_) J 2 n (k±a) kj_dk ± , (A3) 

Jo 4 J 

/>oo />oo 

* n (fc||,a) = / V n+ i (fe|| , s, a) sds = / V n {k\\ , a, s) sds (A4) 
Jo Jo 



( 27 r) 4 



U 2 (k h k±) Jl{k ± a) kldk±. 



4 Jo 

Using the recurrent relations between the Bessel functions we obtain 

2n d 



*„_i (fcii, a) - * n+ i (fcii, a) = — — (An, a) . (A5) 
Thus the integrated energy transfer is expressed only by the functions $„ 

°° ( g s 2 d 
x ^ n 2 I 3$„ (fcii, a) +fcii-r^$„ (fcii, a) H — $„ (k\\,a) 



a da 
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As an important particular case consider the regularized and screened interaction potential Eq. ([52)) . For this potential 
Eq. flMJ) for the function yields [H 

2 11 

*n (k , a) = A _ 2 _ d _ 2 [In (ko) if„ (ko) - /„ (xa) if„ (xa)] - ^ #„ (ko) - ~^i Bn (x a ) > (A7) 

where k, x and d have been introduced in Sec. IIVI and 

H n (0 = ^ [/« (0 x n (0] = i [i' n (0 x„ (0 + In (0 < (0] ■ (A8) 

Note that i?„ (£) = T„ (£,£), where T n (x,y) is defined in Eq. (|55p. From Eq. (|A7|) one can derive the function $„ 
for screened and Coulomb potentials. For the screened but unregularized potential (i.e. A — > 0) all terms in Eq. (|A7|) 
containing x vanish and the function $„ becomes 

$ n (fe|| , a) = - TnHn ( Ka ) • ( A9) 

In a limit A — > oo, i.e. in the case of unscreened Coulomb potential in Eq. (|A9j) the variable k is replaced by fey I. 



APPENDIX B: INTEGRATED ENERGY TRANSFER FOR AN INFINITELY STRONG MAGNETIC 

FIELD 



Consider the integrated energy transfer for an arbitrary ion motion and for infinitely strong magnetic field. The 
integration of Eq. (|69|) with respect to s involves two integrals of the functions T 2 2 (s) and Xbi (s) T03 (s) which can 
be evaluated employing the relation (|A1[) . In a general case with g ^ we obtain 



poo poo (27T"1^ Z" 00 

T (?,*)=/ T 1 2 2 (q, S ) S ds= T 01 (q,s)T Q3 (q,s)sd S = ^-L JJ 2 (q, k ± ) k 3 ± dk X . (Bl) 
Jo Jo 4 Jo 



Obviously the integrations in Eq. (|Bip require that the interaction potential must decay faster than r 1 at large 
distances and must increase slower than at small ones. In particular, for regularized screened potential from 
Eq. jBl|) we find 

2g 2 A 2 + x 2 + l , g 2 A 2 + x 2 
2 (h a — 1) g^A z + 1 



with x = X/d = 1 + X/X which at q — is simplified to 



Uq(x) — T (0, h) = —TT^—r ln x — 1. (B3) 
?r — 1 



Thus the s-integrated energy transfer in the presence of an infinitely strong magnetic field reads 

2n AAif %sds = 2lTZ2 ^ Uo (x) fe,, - t, 2 - . (B4) 

Jo 

The quantity W (x) in Eq. (IB4j) can be treated as a modified Coulomb logarithm. 
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